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The issue of a self-consistent solution of Maxwell-Einstein equations achieves a very 
simple form when all quantum effects are neglected but a weak vacuum polarization 
due to an external magnetic field is taken into account. From a semi-classical 
point of view this means to deal with an appropriate limit of the one-loop effective 
Lagrangian for electrodynamics. When the corresponding stress-energy tensor is 
considered as a source of the gravitational field a surprisingly bouncing behavior 
is obtained. The present toy model leads to important new features which should 
have taken place in the early universe. 

1 Introduction 

Classical electrodynamics provides results in very good agreement with cosmic ob- 
servations from now up to the epoch near the end of the radiation era. During 
the radiation dominated epoch, however, it is commonly adopted some inflationary 
mechanism, whose origin relies on quantum theory: inflatons, strings, membranes, 
etc. With this tool one describes the early universe matter content as an ionized 
plasma together with a large scale magnetic field,EI whose origin may rest on cos- 
mic strings ,□ pseudo Goldstonc bosonsa or other sources. At the GUT scale this, 
magnetic field may reach very high intensities, H ~ 10 40 Gauss or even stronger .□ 
Fields of this strength by far exceed the limit H c pa 4.4 • 10 13 Gauss beyond which 
QED vacuum polarization must be considered. Onti-loop effective potentials were 
thus calculated for a background magnetic fieldJ§LrQ the influence of them on the 
primordial magnetohydrodynamics being recently proposedB 

Nevertheless, no attempt to analyze the back-reaction of such corrections on 
the gravitational field itself was devised. In this vein, we will here limit our con- 
siderations to a toy model which describes the weak field limit of the one-loop zero 
temperature effective Lagrangian of QED driven by an external magnetic field.Q 
The simplicity of the model allows its full solvability, yielding for a Friedmann- 
Robertson- Walker (FRW) spacetime a bouncing solution. A natural upper bound 
for the magnetic field then arises, thus providing its mathematical consistence (i.e., 
the existence of a non-singular magnetic field defined throughout the whole history 
of the universe) . The rather distinct behavior of the above solution as compared 
with its classical counterpart (for which an initial singularity is unavoidable) sug- 
gests that the semi-classical treatment of the matter content of the actual universe 
could provide a very powerful tool for dealing with a cosmic singularity (and all the 
difficulties of standard cosmology from it derived, as the horizon problem). Other 
gauge interactions could similarly be considered, the cosmological relevance of them 
occurring as well at energy scales for which the Planck mass can still be neglected. 

In section we trace the origin of cosmic singularities in the most simple case 
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of a FRW universe in the radiation era. It also presents the required spatial mean 
value algorithm to make such a model consistent with its own isotropy. Section |^ 
applies this same procedure to the weak field limit of QED effective Lagrangian, 
whose dynamics is obtained and solved in section 3.1 . Finally, section [| shows 



that classical ultra-relativistic matter fields (being the extreme limit of fast moving 
ions of the primordial plasma) cannot modify the regularity of the above referred 
solution. 



2 Einstein-Maxwell singular universes 



Classical Maxwell electrodynamics leads to singular universes. In a FRW framework 
this is a direct consequence of the singularity theorems Jj which follows in this case 
from the exam of the energy conservation law and the Raychaudhuri equation for 
the Hubble expansion parameter 9. Let us set for the line element the form 

ds 2 =c 2 dt 2 - A 2 (t)da 2 . (1) 

The 3-dimensional surface of homogeneity is orthogonal to a fundamental class of 
observers endowed with four- velocity vector field = . In terms of the scale- 
factor A(t), the expansion parameter is defined as 

9 = 3 A/ A, (2) 

where a 'dot' means partial derivative with respect to time (i.e., Lie derivative with 
respect to the velocity v). 

For a perfect fluid with energy density p and pressure p the energy conservation 
law and the Raychaudhuri equation assume respectively the form 

p+{p + P )9 = 1 (3) 



e + \o 2 = -^{p + i P ), (4) 

in which k = 8irG/c 4 is the Einstein gravitational constant. Equations (0) and (|J) 
do admit a first integral 

9 2 = 3np. (5) 

Equations (|J) and (||), which also imply (||) as well, can be identified with the 
timelike and the trace components of Einstein equations. 

Since the natural spatial sections of FRW geometry are isotropic, electromag- 
netic fields can generate such universe only after a suitable spatial average be per- 
formed lii The standard procedure^ is just to setf] for the electric field Ei and the 
magnetic field Hi the following mean values: 

< Ei > = 0, (6) 

<Hi> = 0, (7) 



a We will restrict ourselves throughout to the exam of the Euclidean section case. 

6 We make use of Gaussian Cartesian coordinates. Latin indices run in the spatial range (x, y, z), 

while Greek indices run in the spacetime range (t, x, y, z). 
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< EiEj > = -^E 2 9ij , (8) 

< HiHj > = -~H 2 9ij , (9) 

< Ei Hj > = 0. (10) 

Here E 2 and H 2 are both nonnegative functions of time^J, and we denote by angular 
brackets the volume spatial average {e.g., < X > represents the volume average of 
the arbitrary quantity X) for a given instant of time t, 

< X >= lim — / Xd 3 x\ (11) 
v^v a V J 

where V = J d 3 x l with x % € a being spatial coordinates, and V stands for the time 
dependent volume of the whole space (which is finite for the closed section). 

The canonical stress-energy tensor associated with Maxwell Lagrangian is given 
by0 

r M „ = F,af% + ^V> (12) 
in which F = Fp,, F» v = 2(H 2 - E 2 ). Equations ©-© imply 

< F lta F a v ^(E 2 +H 2 ) Vf ,v„ + ^(E 2 -2H 2 )g^. (13) 

Using this result into the expression ( [1^ ) of the stress-energy tensor, it follows that 
its average value < T^ u > reduces to a perfect fluid configuration 

< >= (p-y +Pj) Vf^Vv -p-yguv, (14) 

with energy density 

P-r = l(E 2 + H 2 ), (15) 

and pressure 

Vi = \pi- (16) 

The fact that both the energy density dig) and the pressure (|l6| ) are nonnegative 
for all time immediately yields the singular nature of classical FRW universes, as 
can be seen from Raychaudhuri equation (ji|). In more precise words, Einstein field 
equations for the above energy-momentum configuration yield for the scale-factor 
the typical behavior 

A(t) - t 1/2 . (17) 

All this is standard and well-known. However, near the maximum condensation 
era, classical Maxwell equations do not provide a correct description of electrody- 
namics. Instead, one should consider its quantum corrections. 



c They are not scalars, however, but depend on the set of coordinates, as far as expression (jllj) is 

not a tensor definition but if X is a scalar. 

d We use Heaviside non-rationalized units throughout. 
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3 Quantum corrections of QED at the radiation era 



The effective action for electrodynamics due tQ.one-loop quantum corrections was 
originally calculated by Heisenberg and Eulerfl The development of this work to 
a gauge-invariant formulation is due to Schwingerfl We present here only the first 
order calculation for the effective Lagrangian density 

L = ~\f + ^ + ^G\ (18) 

in which G = F* V F^ = \v a ^ v F a ^F^ = -4E ■ H, with F* v = ^^pF^. By 
flpvafi we denote the Levi-Civita skew tensor, and 

2 / h \ 3 l 

fi = — a 2 [ w 1.67- 1CT 31 cm 3 /erg, (19) 

45 \m e c J m e c z 

where a w 1/137 is the fine-structure constant. 

Note that the homogeneous Lagrangian (|l8|) requires some kind of spatial av- 
erage over large scales, as given in (^)-(Eo|). If one intends to make similar calcu- 
lations on smallep-acales, then either the more complex non homogeneous effective 
QED-LagpangianllHl should be used or else some additional magnetohydrodynamical 
effectt^EJ should be devised in order to achieve correlationlla at the desired scale. 

Treating such quantum correction as a mere effective contribution to classical 
field theory, the corresponding modified stress-energy tensor readal3 

T IMl/ = -4L F F tl a F a „ + (L G G-L)g IMl/ , (20) 

in which L p represents the partial derivative of the Lagrangian L with respect to 
the invariant F, and similarly for the invariant G. 

Since we are interested mainly in the analysis of the behavior of this system in 
the early universe, for which-thiipactual ponderable matter should be identified with 
a primordial ionized plasmaO'HL3, we are led to limit our considerations hprg t/,i the 
case in which only the average of the squared magnetic field H 2 survives^JijO'EiEj 
This is formally equivalent to put E 2 — in (||). Let us see what the consequences 
of this result are. 

3.1 Equation of Motion and Energy Distribution 

Since the average procedure is independent from the equations of motion of the 
electromagnetic field we can use the above formulae (§)-© with E 2 = to arrive 
at a similar expression as (^) for the averaged stress-energy tensor, again identified 
with a perfect fluid with energy density p 7 and pressure p 7 , which are given by 

p T = iff 2 (l-2 Mj ff 2 ), (21) 
p, ( = iff 2 (1-10 vH 2 ). (22) 



e This is strictly true for a viscosity free ionized plasma. When plasma viscosities are considered 
the resulting mean squared electric field may be non zero, but it would still be much smaller than 
its magnetic counterpart. 
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Since the averaged effective stress-energy tensor is not trace-free, the equation of 
state p-y = p 7 (p 7 ) is no longer_given by the Maxwell prescription (|l(]), but presents 
instead a new, quintessential, Ej term which is proportional to the constant p as 

V~i = \ > P 1 - \v H ^- ( 23 ) 
Equation ( p3|) can also be written in the form 

p y = ~ pj - i [(1 - 2p Pl ) - y/1 - Ap Pl \ 

«^(1-WtK> (24) 

where the classical limit /ip 7 <C 1 was applied to obtain the last equality in (p4|). 

We note that, as p is a positive constant, one could envisage the possibility 
that both the energy density and the pressure could become negative. We shall 
see below that this does not occur for the energy density, but it is precisely the 
case for the pressure. Specifically, there exists an inflationary epoch in the model 
of the universe presented here, for which the inertial energy condition (p + p > 0) 
is violated; the gravitational energy condition (p + 3p > 0) is violated as well. 

Equations (||) and (||) do encompass the whole dynamics. Indeed, from (pl|)- 
( p2| ) we have: 

H 2 (l-4pH 2 ) (f + 2^) -0, (25) 



and 



^ = ^H 2 (l-2pH 2 ). (26) 



Equation (^5|) furnishes^] 



#=5: (27) 

where H Q is a nonnegative constant. With this result, (|2^) turns out to be an 
ordinary first order nonlinear differential equation for the scale-factor 



6A 2 V A 4 



whose solution is| 



,!~ //„-\/|( KC 2 i 2 + 3/x). (29) 



Equation (U7Y) then yields the magnetic field as 



' We shall not consider here the H 2 = constant case, as it does not fit the present behavior of the 
actual universe. 

9 We omit here an integration constant, which represents the origin of the time marks, by arbitrarily 
setting A\ t=0 = A min . 
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whose maximum value 



H„ 



1 



6.14- 10 15 Gauss 



(31) 



is much smaller than the typical values which occurr at the GUT scalefl 

Let us make some comments on this solution. First of all we realize it is a 
bouncing solution^], as displayed in Figure |[ whose 
given by 



A 



(32) 



The actual value of A m i n then depends on the constant H a , which turns out to be 
the unique free parameter of the present model. 




^ure 1. Non-singular behavior of the scale-factor A(t). A m i n is given from @ and t c from 



The energy density p 7 attains its maximum value 

Pmax = 77— ~ 3.75 • 10 27 erg/cm 3 
16 fj, 

at the instant t = t r , where 



(33) 



t r = - \ — ~ 1.64- sec. 



h Fo 
e.g 



C \ K 

th bouncing FRW 



(34) 

solutions see, 
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For lower values of t the energy density decreases, vanishing at t = 0. At the 
same time the pressure becomes highly negative, as displayed in Figure |^. Only for 
times comparable to (or smaller than) t c , which lies beynnd the observational lower 
limit t > t n ~ 10° sec provided by the nucleosynthesis E3, the quantum effects are 
important. Indeed, solution (J2 - 
case at the limit of large times 



fits the standard expression (h7n of the Maxwell 




Figure 2. Time dependence of the electromagnetic energy density p 7 and pressure p 7 . 
given from (|33|) and t c from (p4h. 



The corresponding maximum of the temperature is given by 

1/4 



c 

8/1 (7 



(35) 



where a is the Steffan-Boltzmann constant. In energy units the chosen value of \x 
(which describes virtual pairs of electrons) yields 



12.2 Me V. 



(36) 



This result is just at the upper limit for which QELD. admits a perturbative expansion 
on the temperature-dependent coupling constantcll e, T^J 3 ~ 10 MeV. We do not 
care of this point, since the vacuum polarization process we are interested in does 
not require that the virtual pairs e + e~ be produced in thermal equilibrium wijth the 
electromagnetic field which generates them. Indeed, standard calculations!^ often 
suppose, as we did, that these virtual pairs are created mostly at zero temperature. 

Temperatures of the order of ([36]) can also induce virtual pairs of heavier parti- 
cles, as the up quark. Since the right-hand side of ([35|) increases linearly with the 
virtual mass/charge ratio this would in fact lead to a cascade effect which would, 
presumably end at the heaviest relevant particle: the top quark, whose mass isf^l 
n%t/m e ~ 3.33 • 10 5 . If one considers this additional contribution the maximum 
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thermal energy ( pq ) would then be shifted to about kBT max w 2.70 TeV, occurring 
near t c w 3.33 • 10" 13 sec, while iJ maa; w 3.02 ■ 10 22 Gauss. 

4 Matter fields at the radiation era 

Let us make another comment concerning the influence of the presence of other 
matter fields in the universe. Beside photons there are plenty of other particles, 
and physics of the early universe deals with various sort of fields. In the standard 
framework they are treated in terms of a fluid with energy density p v , which satisfies 
an ultrarelativistic equation of state p v = p u /3. Adding the contribution of this 
kind of matter to the averaged stress-energy tensor < T^ v > of the electromagnetic 
field given by (|l4|) and ( J2l|) and ( p2|) it follows, as usual, that p v is proportional to 
the inverse of the fourth power of the scale-factor 

Pv = pV>A~\ (37) 

This result allows us to treat such extra matter fields as a mere reparametrization 
of the constants H a and fx into H a and /t, given by 

H 2 = H 2 +2pi°\ (38) 
The net effect of this is just to re-scale the value of A m i n as 

1/4 

,„) j Amin- (40) 



\H* + 2pi 

Therefore, it turns out that the phenomenon of reversing the sign of the acceleration 
parameter 6 due to the high negative pressure of the magnetic field is not essentially 
modified by the presence of the ultrarelativistic gas. Only a fluid possessing energy 
density p exo tic which scales as Pexotic = Pexotic ^ n with n > 8 could bepable to 
modify the above result. However, this seems to be a very unrealistic case cl. 

5 Conclusions 

Heisenberg and EuleiS have calculated the effective Lagrangian density to deal 
with the nonlinear electro dynamic effects induced by virtual electron-positron pairs. 
This is valid for frequencies v <C m e c 2 jh. One should wonder if the use of such a 
correction in the framework of Einstein general relativity, as we did before, belongs 
to the above range of applicability. Let us write the on-shell Heisenberg-Euler 
effective Lagrangian ( ]lS| ) in the form 

L = -V(l-2 M i? 2 ). (41) 



The limit of validity of this expression consists in the range 

l-2pH 2 >0. (42) 
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In the history of the universe described above the spatially averaged magnetic field 
strength H 2 is globally regular, and is bounded from above at the precise value for 
which the equality in JI2]) holds, as seen from equation (^l]) and Figure |[ Such 
equality is, however, an extreme limit for the weak field expansion which is assumed 
to hold from the very beginning. The involved series is convergent for all times but 
for the instant of maximum condensation, t = 0. Even worser, the dropped terms in 
the effective Lagrangian ( |l8| ) are not negligible at t = t c also. Therefore, the present 
theory must be regarded as a toy model on cosmology, although being a regular and 
consistent solution of Maxwell-Einstein equations. Moreover, this model explicitly 
provides an example of the cosmological relevance of a semi-classical description of 
matter at the early universe. 

From what we have seen, the cosmic singularity of standard cosmology is linked 
with the strictly classical framework in which both the gravitational field and the 
matter content of the universe are treated.cJ While the temperature increases the 
quantum effects become important. The present paper thus analyzes the weak 
field limit of the one-loop zero temperature electrodynamics from a semi-classical 
approachQ in the realm of a spatially homogeneous and isotropic cosmology driven 
by an external magnetic field.tj Instead of being conclusive, however, the regular 
solution (|3^) suggests that more accurate descriptions of matter fields at the early 
universe may provide as well a globally consistent solution. The scenario devised 
here deserves therefore further investigation. 
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